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Abstract. We introduce analytical expressions for a pseudo fully analytical elliptical projected Navarro, Frenk & 
White (NFW) mass prohle to be used in lensing equations. We propose a formalism that incorporates the ellipticity 
into the expression for the lens potential, producing a pseudo-elliptical mass distribution. This approach can be 
implemented to any circular mass profile for which the projected mass profile E(r) and the deflection angle profile 
a(r) both have analytical expressions; however the potential does not necessarily need to take an analytical form. 

We apply this new formalism to the NFW mass distribution and study how well this pseudo-elliptical NFW model 
describes an elliptical mass distribution. We conclude that the pseudo-elliptical NFW model is a good description 
of elliptical mass distributions provided that the ellipticity of the projected mass distribution is < 0.4, although 
with a slightly boxy distribution. 
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1. Introduction 

Cosmological N-bod y simulations of cluster formation 
( [Navarro et al. 1997| ) indicate the existence of a universal 
density profile for dark matter halos, independent of their 
mass, power spectrum of initial fluctuations or cosmologi¬ 
cal parameters. For this so-called NFW profile, the density 
increases near the centre with a shallower slope than an 
isothermal profile, while it steepens gradually outward and 
becomes steeper than isothermal far from the centre. Its 
analytic expression is given by 


p(r) = 


(r/rs)(l + r/rsf 


( 1 ) 


where pc is a characteristic density and a scale radius. 


Recent higher-resolution simulations (e.g. Moore et al. 
1995 ; Ghigna et al. 2000| ) advocate a steeper central cusp 
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of p oc Attempts to constrain the inner slope of 

the density profile with high resolution observations of lu¬ 
minosity profiles ( [Faber et al. 199'? ) seems to confirm a 
central cusp {p cx r“^), rather than a core radius for mas¬ 
sive galaxies. On larger scales. Smith et al. (2001) used 
gravitational lensing to constrain the density profile of 
A 383, a massive galaxy cluster at z = 0.19, finding a log¬ 
arithmic slope of ~ —1.3. Robust interpretation of these 
observational results is complicated by several factors, in¬ 
cluding the absence of baryons from high resolution nu¬ 


merical simulations, systematic uncertainties in the lens 
models arising from parametrisation of the mass distribu¬ 
tion, and the need to use elliptical mass distributions to 
fit observed multiple image systems. 

Gravitational lensing is an ideal tool to constrain the 
radial structure of collapsed halos such as galaxies and 
clusters of galaxies ( Smith et al. 2001 ). However, lensing is 
only sensitive to the projected mass distribution, and ellip¬ 
tical mass distributions are needed to match the multiple 
images observed in both galaxy and cluster lens systems 
(Kneib 2001). In response to the debate regarding the in¬ 
ner slope of the density profile, Munoz et al. (2001) intro¬ 
duced a general set of ellipsoidal lens models with p oc r~^ 
as r —!■ 0 and p oc r“" at large radius. However, as there 
are no general analytic expressions for cusped ellipsoidal 
mass models, the deflections and magnifications are calcu¬ 
lated numerically. They applied their model to the gravi¬ 
tational lens ARM 08279-1-5255 and found a very shallow 
cusp (7 < 0.4). In contrast, for B 1933-1-503, they found 
that a steep density cusp (1.6 < 7 < 2 .0) is favoured . 
To avoid expensive numerical integration, Barkana (1998 ) 
suggested an alternative solution. For a softened power- 
law elliptical mass distribution, it is possible to approx¬ 
imate the integrand so that the integration can be done 
analytically. Therefore, for this flat core model, the deflec¬ 
tion can be then calculated to high accuracy. 

In this paper we propose a new way to introduce el¬ 
lipticity in lensing model in a fully analytical way, and 
we discuss in detail the recipe and limit of the model for 
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the NEW mass distribution. In Sect. 2, we briefly discuss 
spherical NEW lens models. Then we present, in Sect. 3, a 
general pseudo-elliptical formalism that incorporates the 
ellipticity in the expression of the lens potential if this is 
known, or anyway of the deflection angle. In Sect. 4, we 
apply this formalism to the NEW profile and study the de¬ 
parture of this model from an elliptical NEW mass model. 
Einally, in Sect. 5 we discuss prospects for the application 
of this new formalism. 


2. Spherical NFW Lensing Model 

We first recall the expressions for the spherical NEW den¬ 


sity profile (e.g. Bartelmann 1996; Wright & Brainerd 
200 C), this will also allow us to define all the lensing quan¬ 
tities used hereafter. 

In the thin lens approximation, we define z as the opti¬ 
cal axis and z) as the three-dimensional Newtonian 
gravitational potential - where r = \/B? . The re¬ 

duced two-dimensional lens potential in the plane of the 
sky is given by (Schneider et al. 1992|): 


ip{e) = 


D-ls 


+ 00 


Doi^Dos 


^{Doi^O^z) dz 


( 2 ) 


where 0 = { 01 , 62 ) is the angular position in the image 
plane. 

The deflection angle a between the image and the 
source, the convergence k and the shear 7 are then simply: 


a{6) = VeifiO) 

= 2 [-m 


iHo) = \\i{o)r 


dOV 

1 

4 


del) \d6id62j ■ 
( 3 ) 


Eor convenience we introduce the dimensionless radial co¬ 
ordinates X = {x\,X 2 ) = R/ts = 6/9s where Os = 
Xs/DoIj- In the case of an axially symmetric lens, the 
relations become simpler, as the position vector can be 
replaced by its norm. The surface mass density then be¬ 
comes 


E(a:) 


+ CXD 

J p{rs X, z)dz = 2 pcrsF{x) 


( 4 ) 


and the mean surface density inside the dimensionless ra¬ 
dius X is 


S(a:) = —- / 2TTx'S{x)dx = ipcV, 


9{x) 




( 5 ) 


with 




1 1 , , 
arcch— {x < 1 ) 

X 


g{x) = < 


1 


•\/l — X 


1 + ln- 

In- -I- —= 
2 


{x = 1) 
arccos— (x > 1 ) 

X 


The lensing functions a., k and 7 also have simple ex¬ 
pressions (Miralda-Escude 1991) 


E(x) 


0 


= 0 ^^^ = 4ks —5r(x)e:i 

Sent x‘‘ 

k{x) = =2ksF{x) 

_ ^crit 

-fix) = - — =2k 


( 6 ) 


Sprit 






with Ks = Noting V^aix) = {d^^^a, d^^a) and 

(j) = arctan(x 2 /xi), we obtain some useful relations for the 
following that hold for any circular mass distribution: 


k{x) = 

7(x) = 

dxiCtjx) _ 

COS(j) 

k{x) + 7 (x) = 


dx 2 ^\X) 

sinc^ 

a{x) 

OsX 


' a{x) ^ dx^ajx )' 
X cos (f) 

'a{x) _ dxiajx )' 

coscj) 


( 7 ) 


By integrati ng the deflection angle , we find the lens 
potential (p{x) ( Meneghetti et al. 2002| ): 


(p{x) = 2 k,s 01 h{x) 


( 8 ) 


where 


h{x) = 


In^ — — arcch^ 
2 
2 X 


(x < 1) 


In — -I- arccos"^ — (x > 1) 

2 X 


( 9 ) 


— C30 


with 


1 


F{x) = < 


x2 - 1 
1 

3 

1 


x^ — 1 




arcch- 



1 1 
, arccos — 
X 


(x < 1) 
(x = 1) 
(x > 1) 


The velocity dispersion a{r) of this potential, com¬ 
puted with the Jeans equation for an isotropic velocity 
distribution, gives an unrealistic central velocity disper¬ 
sion cr(0) = 0. In order to compare the pseudo-elliptical 
NEW potential with other potentials, we define a scal¬ 
ing parameter Vc (characteristic velocity) in terms of the 
parameters of the NEW profile as follows: 

= ^Gr^Pc ( 10 ) 
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Fig. 1. System of 5 multiple images generated by a pseudo-elliptical NFW cluster at = 0.3 with the following lens 
parameters: Vc = 2000 km/s, Os = 31.3"(rs = 150 kpc) and different values of e. From left to right : e = 0.1, 0.2, 0.3. 
Solid lines are the critical and caustic lines for a source redshift zsi = 1. Dashed lines are the iso-contours of the 

dimensionless projected density —. Units are given in arcseconds. 

^PcTs 


Using the value of the critical density for closure of the 
Universe pcrit = SHq/SttG, we find 


Pcrit 


■^0 ' s 


= 1.8 Wh 


,3 u-2 


150 kpc 


-2 


2000 km s 


( |Golse 


We showed 

2000 km s“^ corresponds to a velocity dispersion 


1200 km s ^ for a Hjorth fc Kneib (2002| ) model. 


— y/die 
X2e = y/ a^X2 

Xe = \J x\^ + X\^ = x/aieXf + a2eX 
= arctan (a: 2 e/xie) 


a.s{x) = 


^ = a(x<;)^^cos(;ie 

(JX\ 

\ ^ a(a:<;)v^sin(/)j 


( 11 ) 


For instance, Meneghetti et al. 
NFW elliptical model: 

aie = 1 - e 

a2e = 1/(1 - e) 


( 2002 | ) use for their 


(14) 


et al. 2002 ) that a value Vc = 


3. Elliptical Potential and Deflection-Angle Model 

We will here introduce an ellipticity e in the circular lens 
potential (p{9). Moreover, we assume that the radial profile 
can be scaled by a scale radius Og, thus making possible 
to define x as x = O/Og (one can always set = 1 if the 
radial profile is scale free). We introduce the ellipticity in 
the expression of the lens potential by substituting x by 
Xg, using the following elliptical coordinate system: 


for an ellipticity along the xi axis. This choice has the 
advantage to stick to a ’standard’ definition: e = 1 — b/a 
- where a and b are respectively the semi-major and semi¬ 
minor axis of the elliptical shape. 

However this choice of au and a 2 e does not yield sim- 
ple expressions for len sing quantities e.g. k and 7 (see 
[Meneghetti et al. 2002 )). Nevertheless, we will now show 
that it is possible to derive simple analytic expressions of 
K and 7 for a particular choice of au and 02 ^. 

At this point, our proposed method can be consid¬ 
ered twofold, i) Either the circular lens potential (p and 
the 2D surface mass density E both have analytic expres- 


( 12 ) 


where au and a 2 e are two parameters defining the ellip¬ 
ticity. 

Furthermore, from the elliptical lens potential ipe (a:) = 
(p{xg), we can then compute the elliptical deviation angle: 


(13) 


sions. We can then introduce the elliptical formalism (12) 
in the lensing potential ip and derive the elliptical deflec¬ 
tion angle cy.g(x) (Eq. ([l3[)). ii) Or, there is no analytic 
expression for the circular lens potential (indeed, in many 
cases the circular lens potential has not a simple analyti¬ 
cal expression). In this case, we need analytic expressions 
for both the circular deviation angle a and the 2D sur¬ 
face mass density E. The elliptical formalism ( |l2[ ) is then 
introduced in the expression of the deflection angle as in 
Eq. dil)- The way the deviation angle is defined ensures 
that a.{x) derives from a lens potential ipe{x) = ip{xg), 
even if there is no analytical expression for ip{x). 

Thus, in the following, we will refer to this method 
as the elliptical deflection angle model, whether the lens 
potential is analytically known or not. To be able to sim¬ 
ply derive the convergence and the shear, we choose the 
following elliptical parameters: 


We note that these expressions holds for any definition 
of Ole and a 2 e- 


Ole = 1 - e 
026 = 1 -I- e 


(15) 
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For small e, it gives the same ellipticity along the xi axis 
as the one given by the parameters defined in Eqs (|l^. 
More generally, if we denote by the ellipticity of the 
lens potential contours - taken as 1 — 6/a we have: 


— 1 



(16) 


independently of the profile. This means there is a direct 
relation between the ’standard’ ellipticity and the deflec¬ 
tion angle ellipticity. 

However, for this particular choice of e we can derive 
easily - using Eqs ([^ - the corresponding convergence 
K^(x) induced by Eq. dil): 


Ke(x) 


1 




dxl 


dx\ 

K{Xf) + ecos2(j)^j{x^). 


e / d^p{x^) _ 


9x1, ) 


(17) 


Similarly, the shear j,{x) can be written as: 




e 





40f y \ dx\ dx^ ) \ dxidx 2 ) J 

= 

/ d‘^ip{x,) _ d‘^p{x,) \ / d‘^p{x,) d‘^p{x,) \ 

I dxl, 9x1, ) V 9x1, 9x1, ) 

j / d^pjx,) d^(p(x,) \ ^ 9^ip{x,) \ ^ 1 

(V 9x1, 9x1, J \dxi,dx 2 ,) j 

(18) 


which, using Eqs (^, can be simplified as: 


^l{x) = "f^{x,) + 2 ecos 24 >,"f{x,)K{x,) 

+ e^{K'^{x,) - cos^ 2(l),j'^{x,)). (19) 

Finally, the projected mass density Tt,{x) is simply 
determined from equations ( 0 ) and (||): 

T„{x) = S(a;£) -|- €COs2(j),{T,{x,) — T,{x,)). (20) 


4. Application of the Model to NFW Halos 

Now, we apply the elliptical deflection angle model devel¬ 
oped in Sect. ^ to the NFW profile (p. In that case, the 
lens potential (i and the 2D projected mass profile (i 
are known analytically. 

An illustration of some lensed images using our new 
formalism applied to the NFW profile is shown in Fig. |^. 
The caustic associated with the tangential critical line has 
the usual diamond shape and is not reduced to a central 
point as in the spherical NFW case. This of course makes 
the formation of 5-image configurations with tangential 
images possible. 



e 

Fig. 2. Ellipticity es of the projected density T,, as a 
function of the ellipticity e defined in Eq. 0) with the 
choice (|^ for the NFW profile. We show curves for dif¬ 
ferent values of r/xs (r: ellipse semi diagonal, rgi NFW 
scale radius). 



e 

Fig. 3. Parameter €^2 from Eq. (^) used to fit the pro¬ 
jected density S^; it characterises the deviation from a 
purely elliptical curve. e-s 2 is a function of the ellipticity 
e defined in Eq.(|l^ with the choice ( 0) for the NFW 
profile. Curves are shown for different values of r/xs {x: 
ellipse semi diagonal, r^: NFW scale radius). 


4.1. Expression of the 3D pseudo-elliptical mass 
distribution 

This particular mass distribution has the advantage that 
the 3D pseudo-elliptical NFW mass profile p,{x,X 3 ) can 
also be derived. Indeed using the scaled variables X 3 = 
z/xs and u = x/xs = \/x^ + x^, we can compute from 
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We now investigate the range of e for which this NEW 
mass model is an adequate description of an elliptical un¬ 
derlying mass distribution. We will use two methods to 
quantify the deviation of our model from a purely ellipti¬ 
cal distribution. 

Fig. shows the contours (dashed lines) of the pro¬ 
jected mass density (Eq. @) for e = 0.1, 0.2, 0.3. In 
the more elliptical models, the contours become increas¬ 
ingly boxy/peanut shaped at larger “radius”. In order to 
investigate this boxy behaviour, we must first quantify 
the ellipticity es of the projected mass distribution E,., 
and then relate this to the ellipticity e of the lens model. 
Purely elliptical projected mass density contours would 
have a polar equation of the type 


1 


'(l-eE)cos2</.+ |^ 


(23) 


es 


r oc 


(1 - es) cos^ (j) + 


sin^ (j) 
1 — es 


(24) 


es2 cos(4(^ej,) 


with 


tani^ 

= arctan -- 

^ ' 1 - es 


(25) 


We propose a fit of elliptical-like functions which is a 
deviation from an elliptical model. It is slightly different 
from the function presented by Jedrzejewski (1987) and 
already used by Shaw (19^ ) or Quillen et al. (1997 ). We 
write the polar equation 


Fig. 5. Goodness of fit of the projected density E^ con¬ 
tours by a function of the type (|^). The ellipticity e is 
defined for the NEW profile in Eq. ©) with the choice 
(ITsI). Curves are shown for different values of r/r^ (r: el¬ 
lipse semi diagonal, rg'- NEW scale radius).The goodness 
of fit is computed as in Eq. ( p^ with N = 20. 

deviation of a boxy/peanut function from an ellipse via 
the parameter es 2 - 

Eor a given ellipticity e, introduced in the deflection 
angle, and a given radius r = va^+^, we will fit the 
parameters es and eE 2 in the corresponding surface den¬ 
sity profile. A goodness of fit indicator will allow us to 
check how effective the representation is. The ratio b/a 
gives a first relation. The other one is given by c/a where 
c is such that Ee(c, ctan ^q;,) = Ee(a,0) = Ee(0, 6). This 
means that we adjust the coefficients of the fitting function 
along the first order ellipse diagonal. Eq. (|^) is indeed a 
deviation from an ellipse in this direction. Actually, get¬ 
ting two relations does not lead analytically to es and e^ 2 , 
mainly because of the angle which depends on es. So 
we assume in practice that 1 — ~ 6/a. This approxima¬ 

tion is correct since 

6 VI — es + es2 /, \ / 2 — es 

- = I . = (1 - es) 1 + - 

all \ 2 

es2 


1 — es 


(26) 


If we assume that the fitted contour is roughly an ellipse 
with an ellipticity es, the angular direction (/ab of its diag¬ 
onal is such that tan (pab = 6 /a~l —e^.In this expression, 
a and 6 are defined by Eg (a, 0) = Eg(0, 6 ), i.e. the “pseudo” 
semi major and minor axis. Taking a radial coordinate of 
the for m (p^ permits to quantify the degree of boxiness 
for a n pn elliptical model. Indeed, compared to an ellipse, 
r is then smaller along the axis and larger along the diag¬ 
onals for es 2 > 0, i.e. the distribution is boxier. This kind 
of fit can be generally applied to check quantitatively the 


for eE 2 ^ 1 and 1 — es = 0(1). It is then possible to ex¬ 
press es and eE 2 analytically for given e and r, see Eigs @ 
and (||). 

We note that a given value of e corresponds to a higher 
value of ce (Fig. (§)). e can be considered as the ellipticity 
of the potential e^p for a large range of values (see Eq. (|^ : 
there is less than 10% error for e < 0.25). It is also known 
that the ellipticity of the projected mass density is pro¬ 
portional to and larger than the ellipticity of the poten¬ 


tial in the linear approximation and then flattens (Kneit 


1993). For instance, a singular isothermal ellipse satisfies 


= Ocp, for Cp, < 1 . 

To derive numerically such a relation for the NFW pro¬ 
file, we need to know the range of acceptable and physical 
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Fig. 4. Solid lines: projected density Y,^{xi/rs,X 2 /rs) contour for a pseudo-elliptical NFW profile at r/xs = 10 (r: 
ellipse semi diagonal, rji NFW scale radius). Dotted lines: best fit ellipse. Dashed lines: fitting function with 
computed parameters es and e-£. 2 - From left to right: e = 0.1,0.2,0.3. 


values for e. For all ellipticities and up to r/r^ = 10 H.e. 
r ^ 1.5 Mpc for a galaxy cluster), es 2 < 0.1 (see Fig. (^). 
So the deviation parameter is not too large and the ellip¬ 
tical approximation could be considered as acceptable if 
the goodness of fit for the function is small. To check 
the relevance of this fit, we plot the contour, the first 
order ellipse and the fitting function found for r/vs = 10 
and e = 0.1,0.2,0.3 (Fig. (^). The fit is correct for small 
ellipticities but is not suited for e = 0.3. In particular it 
fails to reproduce the shape along the X 2 axis. 

Quantitatively, we define a goodness of fit in the fol¬ 
lowing way: 


gof 


1 - ri(j)i)\ 

N ^ r((/)i) 


(27) 


TT Z 

for ~ '2 symmetry reasons). For given a and e, 

and r{(pi) are respectively the distances of the pro¬ 
jected density contour and of the corresponding fit func¬ 
tion ( p^ ) from the centre in the direction (j>i. Fig. i) con¬ 
firms that the goodness of fit becomes quite large from 
e ^ 0.25 (the deviation from the proposed function then 
reaches 10%). 

We think that function (|^ can be useful to test de¬ 
viation from ellipticity of a given function in various sets 
of problems. In our case the deviation parameter es 2 is 
rather small but the goodness of fit is only acceptable for 
ellipticities (introduced in the deviation angle) of e < 0.25. 

Alternatively, to simply quantify the degree of boxi- 
ness for this pseudo-elliptical NFW model, we defined the 
characteristic deviation from ellipticity in the following 
way. On Fig. |^, 5r is the distance between a real ellipse 
and a contour along the ellipse diagonal. We plot 5r/r 
versus e for different t/ts ratios in Fig. ^ At all radii, 
and for all e, the model has a positive 5r, i.e. the model 
mass distribution is more boxy than an elliptical distri¬ 
bution. Assuming that the underlying mass distribution 
is elliptical, and aiming to incur an error in r which is 
< 10%, we find that on scales of 1.5 Mpc (i.e. correspond¬ 
ing to r/xs ^ 10 for a galaxy cluster), the pseudo-elliptical 



Fig. 6. Method used to compare a projected density con¬ 
tour and a real ellipse with semi axes a and b. 5x/x char¬ 
acterises this deviation. 


model provides an adequate description of the underlying 
mass distribution for e < 0.25, which translates to a limit 
of es ^ 0.4 on the projected density at x/xg = 1 (see 

Fig.|). 

For models in which the potential - rather than the 
deflection angle - is chosen to have elliptical contours, 
the corresponding density contours acquire the artificial 
feature of a dumbbell shape, and the density can also be¬ 
come negative ( Kassiola &: Kovner 1993 ). Similarly here, 
for large ellipticities or at large radii, we see from Eq.(|^ 
that the projected density can also become negative. 
This occurs closer to the centre along the X 2 axis where 
cos2(^g = — 1. For each value of the ellipticity e, we plot 
in Fig. H the scaled distance b/xg at which Ee(0,6) be¬ 
comes negative. If we decide to have physical (i.e. positive) 
mass density up a scale of 1.5 Mpc (typically b/xg = 10 
for a cluster), we have to restrict ourselves to ellipticities 
smaller than e ^ 0.3 (i.e. es ^ 0.6 at x/xg = 1 from Fig. ||). 
So a relatively broad range of systems can be modelled in 
a physically consistent way. 
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e 

Fig. 7. (5r/r (as defined in Fig. ^ as a function of e. It 
characterises the deviation of the projected density from 
an ellipsoidal model for various r/r^ ratios (r^: NFW scale 
radius). 



Fig. 8. h/rg {b: distance from the centre along the X 2 axis 
at which becomes negative, NFW scale radius) as 
a function of e. 


We want to obtain an explicit, if approximate, relation¬ 
ship between the ellipticity e introduced in the deviation 
angle and the ellipticity of the projected mass density es 
it induces. In the acceptable and physical range [0,0.25] 
for e, we fit a polynomial of the form: 

ce = oi e -I- 02 e^. (28) 

A fit for r = Ts leads to es = 2.27 e — 2.03 with a 
= 3.9 10“”^. More generally, the coefhcients Oi depend 
on X = r/rg. A fit between x = 0 and x = 10 gives 


r oi = 2.12 + 0.179X . . 

I a2 = -1.70 - 0.328 X 

with a = 2.9 10"^. 

In summary, we can say that the deviation angle el¬ 
liptical model can be applied to NFW mass profile up to 


e ^ 0.25. For this range of values, e can be identified with 
the ellipticity of the potential e^, and the ellipticity of the 
projected mass density ce is about twice larger than e. 

5. Conclusion 

We propose a simple new formalism that introduces the 
ellipticity into the lens potential/deflection-angle of a cir¬ 
cular mass model. The method can be applied when the 
lens potential or/and the deviation angle takes an ana¬ 
lytical form. Then for radial mass profiles for which the 
2D surface density E also has an analytical expression, 
this formalism gives analytical expressions of a pseudo¬ 
elliptical mass distribution for the deviation angle, the 
projected mass density, the convergence and shear. 

Whatever the form of the mass distribution, the el¬ 
liptical parameter e is simply expressed as a function of 
the ellipticity of the potential. This is particularly helpful 
in getting some insight on the physical meaning of this 
parameter. 

We have applied this formalism to the NFW profile and 
estimated the range of ellipticity (e < 0.25, or < 0.4) 
for which this model is a good description of elliptical 
mass distributions and thus can be reliably applied to ob¬ 
servational data. To derive these limits, we introduced a 
particular fit for elliptical-like profiles, that can be useful 
in similar cases. 

Our proposed method is particularly useful when it 
is essential to quickly calculate the potential, the deflec¬ 
tion angle and magnification of many images and/or many 
mass clumps. This is particularly important when using 
inverse methods (such as maximum likelihood) to inves¬ 
tigate galaxy-galaxy leasing in the field or in clusters of 
galaxies, or to compute time delays. 
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